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Abstract 

We analyze ill-defined pinch singularities characteristic of out of equilibrium thermal 
field theories. We identify two mechanisms that eliminate pinching even at the single self- 
energy insertion approximation to the propagator: the first is based on the vanishing of 
phase space at the singular point (threshold effect). It is effective in QED with a massive 
electron and a massless photon. In massless QCD, this mechanism fails, but the pinches 
cancel owing to the second mechanism, i.e., owing to the spinor/tensor structure of the 
single self-energy insertion contribution to the propagator. The constraints imposed on 
distribution functions are very reasonable. 
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1 Introduction 



Out of equilibrium thermal field theories have recently attracted much interest. From the 
experimental point of view, various aspects of heavy-ion collisions and the related hot QCD 
plasma are of considerable interest, inparticular the supposedly gluon-dominated stage. 

Contrary to the equilibrium case''O0^ where pinch, collinear, and infrared problems have 
been successfully controlled'-B S B B-*, out of equilibrium theory^ 0' B 0) has suffered from 
them to these days. However, progress has been made in this field, too. 

WeldonlI9 has observed that the out of equilibrium pinch singularity does not cancel; hence 
it spoils analyticity and causality. The problem gets worse with more than one self-energy 
insertions. 

Bedaque has argued that in out of equilibrium theory the time extension should be finite. 
Thus, the time integration limits from — oo to +cx). which are responsible for the appearance 
of pinches, have to be abandoned as unphysicaldl. 

Le Bellac and MabilatHH have found that collinear singularities cancel in scalar theory, 
and in QCD using physical gauges, but not in the case of covariant gauges. Niegawa0 has 
found that the pinch-like term contains a divergent part that cancels collinear singularities in 
the covariant gauge. 

In their discussion of the pinch-like term Le Bellac and Mabilat0 try to avoid the problems 
with pinching singularity by substituting the bare retardedphoton (gluon) propagator with the 
resummed Schwinger-Dyson series calculated in the HTL0. 0. approximation: 

with the spectral function given by pt,l = 2/mAT,L(Q'o+^^, <?)• In the above expressions m is the 
thermal photon (gluon) mass given by = (m^ = g'^T'^{Nc + Nf/2)/6 for the gluon; 

note that they assume a small deviation from equilibrium), and x = qo/q ■ In equilibrium, 
at high temperatures and low momenta, this substitution is necessary in order to obtain the 
results correct to the leading order in g'^. One expects that similar methods work also for a 
narrow class of particle distributions corresponding to "high temperatures" out of equilibrium. 
However, their expression for pinch-like contribution differs from the one following from our 
general expression for the resummed Schwinger-Dyson series ( |2.32D . At medium and large 
photon momenta the elimination of pinching by the use of the HTL approximated propagator 
is no more justified. 

Altherr and Seibert have found that in massive g'^(j)^ theory pinch singularity does not 



occur owing to the kinematical constraint0. This result is restricted to the case of one-loop 
self-energies. 

Altherr has suggested a regularization method in which the propagator is modified by the 
width 7 which is an arbitrary function of momentum to be calculated in a self-consistent way. In 
g^(j)'^ theory, for small deviations from equilibrium, 7 was found to be just the usual equilibrium 
damping rate@. 

This recipe has been jus tified in the resummed Schwinger-Dyson series in various problems 



with pinchingGi' 13' H HI HI . 



Baier, Dirks, and Redlichli^ have calculated the n — p self-energy contribution to the pion 
propagator, regulating pinch contributions by the damping rate. In subsequent papers with 
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SchiffS 

they have calculated the quark propagator within the HTL approximationliS LL21 1121 ; 
in the resummed Schwinger-Dyson series, the pinch is naturally regulated by ItjiLr. 

Carrington, Defu, and Thoma@ have found that no pinch singularities appear in the HTL 
approximation to the resummed photon propagator . 

NiegawaS has introduced the notion of renormalized particle-number density. He has 
found that, in the appropriately redefined calculation scheme, the amplitudes and reaction 
rates are free from pinch singularities. 

By pinching singularity we understand the contour passing between two infinitely close 
poles: 

7 ^7 (1-2) 

where x = — m^. It is controlled by some parameter, e.g., e. For finite e, the expression 
is regular. However, when e tends to zero, the integration path is "pinched" between the two 
poles, and the expression is ill-defined. Integration gives an contribution plus regular terms. 
By performing a simple decomposition of [x ± ie)^^ into PP{l/x) =F iTT6{x), one obtains the 
related ill-defined 5^ expression. 

The following expression, which is similar to ( p..2[ ), corresponds to the resummed Schwinger- 
Dyson series: 

^^7 r-. (1-3) 

{x - J:r{x) + ie) {x - 1:*r{x) -ie)' ^ ^ 



where uj{x) and uj{x) (which appears in ( |1.4| )) are, respectively, proportional to ^{x) and 



where Q,{x), and ^l{x) are the components of the self-energy matrix to be defined in 

Sec. HI. 

rrqi FJTl, pTTl [7J7^ ^T!% 

In expression ( |1.3| ), pinching is absentli^' IfU' IfiJ' Ifrl' If^ if /mS/j(xo) 7^ at a value of Xq 



satisfying Xo — ReT^R^Xo) = 0. 

The expression corresponding to the single self-energy insertion approximation to the prop- 
agator is similar to (|1.3| ): 

dx- "^[^^ . , . (1.4) 

[X + ie)[x — te) 



One can rewrite the integral as 



If it happens that 



dx / 1 1 \uj(x 



2 \x + ie X — iej x 



;i.5) 



lim^^ = K<oo, (1.6) 

then the integral ( |1.5|) decomposes into two pieces that, although possibly divergent, do not 
suffer from pinching. 

There are two cases in which the function uj{x) is even identically zero in the vicinity of 
the a; = point: in thermal equilibrium, because of detailed balance relations; in massive 



^f^^^ theory out of equilibrium, owing to the mass shell condition 113 . The latter mechanism 
also works in out of equilibrium QED if a small photon mass is introduced. However, this 
elimination of pinching can be misleading: the domain of x, where a;(x) = 0, shrinks to a point 
as 0. We shall show that the elimination of pinching also occurs in the rriy = case. 
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In this paper we identify two mechanisms leading to relation (|1.6|) . They are based on the 
observation that in the pinch-like contribution loop particles have to be on mass shell. 

The first mechanism is effective in out of equilibrium QED: in the pinch-like contribution to 
the electron propagator, phase space vanishes linearly as x ^ . In the pinch-like contribution 
to the photon propagator, the domain of integration is shifted to infinity as x ^ 0. For 
distributions disappearing rapidly enough at large energies, the contribution again vanishes 
linearly in the x — limit. This mechanism is also valid in QCD in the cases with massive 
quarks. 

In out of equilibrium massless QCD, phase space does not vanish, but there is an alternative 
mechanism: the spinor/tensor structure in all cases leads to relation ( |1.6|) . 

Also, in out of equilibrium massless QCD, introduction of a small gluon mass does not help. 
In this case, processes like qq ^ g are kinematically allowed, the spinor/tensor structure is 
modified, and VL does not vanish in the x — >■ limit. 

In a few cases, none of the mentioned mechanisms works and one has to sum the Schwinger- 
Dyson series. This is the case of the vr— p loop in the vr self-energy . Even in the limit of zero pion 
mass, uo{x) vanishes only as |x|^''^ and relation ( p.. 61) is not fulfilled. A similar problem appears 
in electroweak interactions involving decays of Z and W bosons, decay of Higgs particles, etc. 
Another important case is massless g'^cfr' theory. In contrast to massless QCD, massless g'^ip^ 
theory contains no spin factor to provide a q^ factor necessary to obtain ( |1.6|) . 

The densities are restricted only mildly: they should be cut off at high energies, at least as 
\ko\~'^~^ 1 in order to obtain a finite total particle density; for nonzero /Cq, they should be finite; 
for ko near zero, they should not diverge more rapidly than |A;o|~^, the electron (positron) 



d have a finite derivative. Further restrictions may come from Slavnov- Taylor 
but they are not crucial for our analysis. 



distribution shou^ 
identities^il 

When necessary we assume that the zero-temperature renormalization has already been 
performed. The "finite temperature" c3 c3 e3 renormalization out of equilibriumED may 
give rise to new problems formally similar to those treated here. Their treatment is beyond the 
scope of this paper. 

The paper is organized as follows. 

In Sec. II we analyze the Schwinger-Dyson equation in the Keldysh representation, solve it 
formally, and identify pinch-like expressions. For one-loop self-energy insertions, we find that 
the Keldysh component {VL{q^)) of the self-energy is responsible for pinches. We further find 
that the nonzero Keldysh component requires loop particles to be on shell. 

In Sec. Ill we analyze functions such as fi, VL, and ItjiLr-, and investigate their threshold 
properties. 

In Sec. IV we show that the electron and photon propagators, calculated in the single self- 
energy insertion approximation, are free from pinching. 

In Sec. V we analyze pinch-like expressions in the g — g, g — g-, and ghost-ghost contributions 
to the gluon propagator, the quark propagator and the ghost propagator in the single self-energy 
insertion approximation. We find that, in all the cases, the spinor/tensor factor F contains a 
factor q^ that is sufficient to eliminate pinching. 

In Sec. VI we briefiy recollect the main results of the paper. 
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2 Propagators and the Schwinger-Dyson equation 

We start @ ^ by defining out of equilibrium tliermal propagators for bosons, in tlie case when 
we can ignore tlie variations of slow variables in Wigner functions 



^=(n" n"), (2.1) 



D21 D22 . 

Dn{k) = DUk) = ^2_J^2re\k^\ + 2vr - m^), (2.2) 

Di2{k) = 2n6{k'^ - m^){cosh.'^ 9Q{ko) + sinh^ eQ{-ko)), (2.3) 

D2i{k) = 2'n5{k'^ - m^){cosh^ eQ{-ko) + sixih^ eQ{ko)). (2.4) 

For particles with additional degrees of freedom, relations ( |2.1| )-( |27^ ) are provided with 
extra factors (Ji + rn) for spin 1/2, {g^^ — (1 — o)k^ky/{k'^ ± 2ieko)) for vector particle, etc., and 
similarly for internal degrees of freedom. To keep the discussion as general as possible, we show 
these factors explicitly only when necessary. The propagator defined by relations ( |2.1| )- (|2^ ) 
satisfies the important condition 

= Z^ii - D12 - D21 + 1^22. (2.5) 
In the case of equilibrium, we have 

sinh^ eiko) = nsiko) = J, . ^ . (2.6) 

exp p\ko\ — 1 

To obtain the corresponding relations for fermions, we only need to make the substitution 

sinh^ e{ko) - sin^ e{ko). (2.7) 
In the case of equilibrium, for fermions we have 

sin^ 9,Ako) = n.Ako) = e,p;5(| ^) + r (^.8) 

Out of equilibrium, UBiko) and npiko) will be some given functions of ko- 
To transform into the Keldysh form, one defines the matrix Q as 

Now 

DA(k) = -D„ + D,, = _ ^ .^^^^ = --DJ(fc) = (2.12) 
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Dxik) = D11 + D22 = 27r5(P -m^)(l + 2sinh2^). (2.13) 
We need Dk expressed through Dr and Da- 

DK = h{K){DR-DA), h{K) = -t{K){l + 2sm\i^e). (2.14) 

Again for fermions, Dk is equal to 

DK{k) =Du + D22 = 27i5{k'^ - m^)il - 2sm^ e). (2.15) 

The proper self-energy 



(2.16) 



satisfies the condition 

= En + S12 + S21 + S22. (2.17) 

It is also transformed into the Keldysh form (in Niemi's paper there is a misprint using 
instead of Q): 

I Q^Q-\ (2.18) 

Sfi = -(Sn + S2i), (2.19) 

S^ = -(Sn + Si2), (2.20) 

fi = Sn + S22. (2.21) 



We also find 



Sa = ^R- (2.22) 



The "cutting rules" (refs.[^, |36|, see also ref. [0 for application of the rules out of equilibrium) 
will convince us that only on-shell loop-particle momenta contribute to IuiLr and VL. 

The calculation of the E matrix gives (propagators S{k) and G{k) in the self-energy matrix 
and in the Schwinger-Dyson equation are also given by (|2.1| ) to (p.l4|) , with the spin indices 
suppressed to keep the discussion as general as possible): 

= -4^?' / -0jiiDR{k){SR{k -q)- SK{k - q)) + {DA{k) - DK{k))SA{k - g)), (2.23) 
= -^\9' J (0i((^^(^) - DK{k))SR{k -q) + DA^SA^k - q) - SK^k - g))), (2.24) 

^ = 4^' / (0^^^^^^^^ + DA{k)){SA{k - g) + SR{k - q)) + DK{k)SK{k - q)). (2.25) 
The Schwinger-Dyson equation 

g = G + iG^g, (2.26) 

can be written in the Keldysh form as 





Qr\ 




Gr\ 






Qk) 


\Ga 


Gk) 





GrErQr 

Ga^aQa Ga^Gr + Gk^rQr + Ga^aQk 



(2.27) 
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By expanding ( 2.27 ), we deduce the contribution from the single self-energy insertion to be 
of the form 

gR^GR + iGr^rGr, gA^GA + tGA^AGA, (2.28) 
which is evidently well defined, and the Keldysh component suspected for pinching: 

gK-GK + iGa^Gr + zGk^rGr + iGa^aGk. (2.29) 

It is easy to obtain a solution® for Qr and Qa using the form ( |2.27| ). One observes that the 
equations for Qr and Qa are simple and the solution is straightforward: 



Gr 



R 



(2.30) 



(2.31) 



To calculate Qk, we can use the solution (|2.3CI| ) for Qr and Qa'- 

Qk = Qa{G^GkG^ + iVL)QR. 

Now we eliminate Gk with the help of (p.l4|) : 

Qk = Qa {h{qo){Gf - Gj,') + tn) Qr. (2.32) 

The first term in ( |2.32| ) is not always zero, but it does not contain pinching singularities! 
The second term in ( p.32|) is potentially ill-defined (or pinch-like). The pinch-like contribution 
appears only in this equation; thus it is the key to the whole problem of pinch singularities. In 
the one-loop approximation, it requires loop particles to be on mass shell. This will be sufficient 
to remove ill-defined expressions in all studied cases. 

Equation (|2.32| ) differs from the one used in Ref. [|12[. Indeed in Ref. |jT2[ Im^^^ is used 
instead of QrQa- 

We start with ( p.29| ). After substituting ( |2.14| ) into ( p.29| ), we obtain the regular term plus 
the pinch-like contribution: 

Qk ~ QKr + Qkp, (2.33) 
QKr = h{qo) {Gr -Ga + iGr-LrGr - iGaT^aGa) , (2.34) 

Qkp = iGa^Gr, n = n- h{qo){^R - s^). (2.35) 

For equilibrium densities, we have S21 = e~^'^°Si2 , and expression ( |2.35| ) vanishes identically. 
This is also true for fermions. 

Expression ( p.35| ) is the only one suspected of pinch singularities at the single self-energy 
insertion level. The function Vt in ( 2.35 ) belongs to the type of functions characterized by the 
fact that both loop particles have to be on mass shell. It is analyzed in detail in Sees. Ill and 
IV (for threshold effect) and in Sec. V (for spin effect). With the help of this analysis we show 
that relation ( p.35|) transforms into 



IKp 



K{q^,w?,qo) 



1 



m? + 2ieqo ~'~ q 



1 



2itqo , 



(2.36) 
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where K{q'^, m?, Qo) is VL/{q^ — im?) multiplied by spinor/tensor factors included in the definition 
of Gr^a- The finiteness of the limit 

hm K{q\m\qo)=K^{qo)<^ (2.37) 

is important for cancellation of pinches. The index =p indicates that the limiting value im? is 
approached from either below or above, and these two values are generally different. To isolate 
the potentially divergent terms, we express the function K{q^, m?, q^) in terms of functions that 
are symmetric {Ki{q'^,m'^,qo)) and antisymmetric {K2{q'^ , 171"^ , qo)) around the value = m^: 

K{q^,m^,qo) = {K^{q\m\qo) + e{q^ - m')K2{q\m\qo)). (2.38) 

These functions are given by 

^i,2(g', m\ qo) = \{K{q\ m\ q,) ± K{2m' - q\ m\ g„)). (2.39) 

Locally (around the value q^ = m^), these functions are related to the limits K±{qo) by 

K,,2{q',m^q,) = ^(ir+(g„) ± i^_(g,)). (2.40) 



As a consequence, the right-hand side of expression (p.36|) behaves locally as 

gKpiq^m^,qo) ^ -UK^{qo) + t{q^-m^)K2{qo)) [- ^Vo- +" 1 o- ) ' 

(2.41) 

and the term proportional to K2 is capable of producing logarithmic singularity. Furthermore, 
we were unable to eliminate pinches related to the double, triple, etc., self-energy insertion 
contributions to the propagator. 



3 Threshold factor 

In this section we analyze the phase space of the loop integral with both loop particles on mass 
shell. Special care is devoted to the behavior of this integral near thresholds. In this analysis 
the densities are constrained only mildly: they are supposed to be finite and smooth, with a 
possible exception at zero energy. We also assume that the total density of particles is finite. 
The expressions are written for all particles being bosons, and spins are not specified; change 
to fermions is elementary. 

To obtain the integrals over the products of Dr a and Sr^Ai we start with a useful relation: 

/+00 \ 
-00 ^'^\k'^ — mj) + iXko€){{k — q^ — nig + iri{ko — qo)^) 

J-oo [k — qY — rrig k^ — mjj 

/+00 
dkof{k, q)e{ko)e{ko - qo)5{e - ml)5{{k - qf - m|) (3.1) 
-00 
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where f{k, g), as a function of ko, is some polynomial of order 0, 1, 2, or 3. 

Relation is obtained as an average of the results obtained by closing the integration 
path through the upper and through the lower semi-plane. 

In the defining relations ( |2.2| ) to ( p.4| ) the particle distribution sinh 9 (and similarly for 
sin^^) appears only in the expression where it is multiplied by 6{k^ — m^). Thus there is 
freedom to replace kg by its on mass shell value e{ko)u!D, where un = (/c^ — m|,)^/^. The physical 
results should not be altered by this replacement. Then we write h{ko), defined in (|2.14|) , as 
—ko/uj£){l + 2 sinh^(a;£))), and similarly for h(ko — g^). This substitution matters when one 
wants to perform integrals in ( p.23| ) to (|2.25|) with the help of With this replacement, 

f{k, q) in (^TTj) does not depend on densities and can be proved as stated above, without 
need to discuss the analytic properties of sinh^ 6 and sin^ 6. Without our replacement, one 
would be immediately stack with the question how these quantities (i.e., sinh^^ and sin^i9), 
supposed to be known, within some uncertainty, along the real axis, behave in the complex 
plane very far from the real axis! 

Similar relations could be obtained for higher powers of Dr^a and Sr^a- For example, for 
the nth power of k"^ — mjj + iXkoe, the real part of the integral will be obtained by substituting 
^(n)^^2 _ ^2^)(_iy") instead of 5{k'^ — mjj). Now we easily calculate -ReS/j as 

_ -^^ r ^4j^^.p 1 - ^Dhpikg) ^ e{ko - qo)S{{k - qf - ml)hs{ko - qo)\ ^ 



2(27r)3 7 \^ ((A;-g)2_m2) (P - mf,) ) ■ 

(3.2) 

F is the factor dependent on spin and internal degrees of freedom. The thermal part of ReT^R 
is given by: 

p V / ^4, 1 5{k^ -ml)sinh]^{ko) <5((fc - g)^ - m|) sinh|(fco - qo) \ 
^'^^'^ = J2^^I'^^^\ ((fc-g)2-m|) + (F - ml) J ^^'^^ 



At equilibrium equation (|3.3|) (after necessary boson— >^fermion conversion) agrees with the 
known results^' . 

Now, starting from ( 2.19| ) to (2.25), we calculate Vt and ImTiR. 



where 



2z/mSn = 2^ / 7^47r25(A;2 - ml)5{{k - qf - ml)Nn{ko, ko - qo)F, (3.4) 



(2vr) 



and 



Nn{ko. ko - qo) = ^i-e{ko{ko - qo)) + (1 + 2 sinh^ OnikoW + 2 sinh^ OsiK - go))), (3.5) 

JmSfi = ^ 1 ^47r25(A;2 - ml)5{{k - qf - ml)Na{ko, ko - qo)F, (3.6) 

Nniko, ko - qo) = (sinh^ 6D{ko)e{ko - qo) + sinh^ 6s{ko - qo)^{~ko) 

+e{-ko)e{ko - qo) - e{ko)&{qo - ko)). (3.7) 
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At equilibrium equations (|3.5| ) to (|3.7| ) agree (after setting F=-l for scalar case) with the 
corresponding equations obtained for boson-boson intermediate state®. 
It is useful to define N^^ko, ko — qo) as 

N^{ko, ko - Qo) = Nn{ko, ko - Qo) - h{qo)Nn{ko, ko - Qo)- (3.8) 

After integrating over 5's, one obtains expressions of the general form 

^ = / ^^^'/'^(^o, ko - qo)F{qo, \q\,ko, \klqk, ...)e(l - z^), (3.9) 
where \k\ = {kl - mj^Y^"^, 

qk = \q\\k\zo, (3.10) 
^ + k^-{q~kf 



2\k\\q\ 

()e(0, 27r) is the angle between vector kx and 

Let us start with the q^ > case. Solution of 9(1 — 2^) gives the integration limits 



(3.11) 



ko 1,2 = ^ {qo{q' + ml- m|) T |g1((g^ - qlr)iq' " Q-tr))'^') , (3-12) 



1^11,2 = (mq' + ml- ml) T qo{{q' - qltM ' q\r)f") , (3-13) 



2q 

or 

"^^"^'^ = i 

q±tr = \mD±ms\. (3.14) 
Assume now that qtr 7^ 0. In this case, at the threshold, the limits shrink to the value 

_ qojql + ml- m|) - _ Iqliql + - ml) 

1^0 tr — „ 2 ' \'^\tr — „ 2 ■ l^O.iOJ 



tr 



Near the threshold, it is convenient to replace the integration variable by (ifco|A;o|/|^| = d\k\. 
Now, for \q^ — qf^\ such small, that the integration limits ko 1,2 are both of the same sign as 

kotr, we have ko = e{ko,tr){k^ + m|,)^/^. 
We define the coefficient Ci by 

ci = ^ j d(j)N{ko tr, ko tr - qo)F{qo, ko tr, \k\tr, qktr, ■■■)■ (3.16) 

Now the expression ( p.9|) can be approximated by 

I ^ ci(|A;|2 - |A^|i) 

^ ci(e(g^ - ql^) + 0( V + dj) '"^^''"''"^^"'-"^^'^' . (3.17) 

Q. 

Relation ( p.l7| ) is the key to further discussion of the threshold effect. 
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We obtain this also for higher dimension (D=6, for example). 

Relation ( |3.17| ) put some limits on the behavior of density functions: they should not tend 
to infinity at any value of go 7^ 0; near Qo = 0, owing to the presence of the factor Qo, they 
should not rise more rapidly than q~^. 

Owing to (|3.14|) and (|3.17] ), the function X(g^,m|,,m|) has the following properties impor- 
tant for cancellation of pinches. 

It vanishes between the thresholds, i.e., the domain (m^) — < < {mj:, + m^)^ is 

forbidden (X = 0). If it happens that the bare mass belongs to this domain, the single self- 
energy insertion will be free of pinching. In this case, multiple (double, triple, etc.) self-energy 
insertions will also be free of pinching. Massive theorylI3 is a good example of this case. 

It is (in principle) different from zero in the allowed domain < {rriD — msY and (m^) + 
msY < q^. In this case, one cannot get rid of pinching. This situation appears in the n — p 
interaction . An exception to this rule are occasional zeros owing to the specific form of 
densities. 

The behavior at the boundaries (i.e., in the allowed region near the threshold) depends on 
the masses m£, and rris and there are a few possibilities. 

If both masses are nonzero and different [0 mo ^ tus 7^ 0), then there are two thresholds 
and X behaves as (g^ — g^j^)^/^ in the allowed region near the threshold g^j^. For m? = g^^, the 
power 1/2 is not large enough to suppress pinching. 

If one of the masses is zero {rriE, 7^ 0, = or niD = 0,ms 7^ 0), then ( p.lTp gives that 
the thresholds are identical (i.e., the forbidden domain shrinks to zero) and one obtains the 
(g^ — mj^Y behavior near mj^. This case (for = mf^) is promising. The elimination of 
pinching in the electron propagator, considered in Sec. IV, is one of important examples. 

If the masses are equal but different from zero {mo = ms 7^ 0), then there are two thresholds 
with different behavior. The function I behaves as (g^ — q^i^Y^"^ in the allowed region near 
the threshold q\^^ = 4m|),and this behavior cannot eliminate pinching in the supposed case 

= 4m|) . 

However, at the other threshold, namely at q'^^^ = 0, the physical region is determined by 
g^ < and the above discussion does not apply. In fact, the integration limits ( |3.12|) or ( p.l3| ) 
are valid, but the region between kg 1 and kg 2 is now excluded from integration. One has to 
integrate over the domain {—00, kg 1) U(^o 25+00). This leads to the limitation in the high- 
energy behavior of the density functions. An important example of such behavior, elimination 
of pinching in the photon propagator (m^), is discussed in Sec. IV. 

If both masses vanish {mo = ms = 0), the thresholds coincide, there is no forbidden region 
and no threshold behavior. The behavior depends on the spin of the particles involved. For 
scalars, the leading term in the expansion of X does not vanish. Pinching is not eliminated. 

The case of vanishing masses {mo = tus = 0) for particles with spin exhibits a peculiar 
behavior. In all studied examples (see Sec.V for details), X behaves as g^ as g^ —>■ 0, which 
promises the elimination of pinching. 



4 Pinch Singularities in QED 



11 



4.1 Pinch Singularities in the Electron Propagator 

In this subsection we apply the results of preceding section to cancel the pinching singularity 
appearing in a single self-energy insertion approximation to the electron propagator. To do so, 
we have to substitute mn = m, ms = 0, sinh^ Oniko) — > —ne{ko), sinh^ Os{ko—qo) n^iko—Qo), 
and h{ko) = —e{ko){l — 2ne{ko)), where Ue and rij are given non-equilibrium distributions of 
electrons and photons in relations (|3.5|) , ( p.7|) , (|3.8| ), and ( p.l4| ). The thresholds are now 
identical 

(4.1) 



and the integration limits are 



ko 1,2 = TTi (io{<f + "^^) =F \q\{{q^ - 



2q^ 



or 



l^li = - m'^))- 



At threshold the limits shrink to the value 

ko tr = Qoy \k\tr = 1^- 



Then, with the help of ( |3.16| ), we define 



If{q\m\qo) 



{q^ — m?) IGvr^lgKg^ — 

(Jl + m)f{qo, \q\,kotr, \k\tr, qhr, •••) W + Hd^b - |^|i)- 
The trace factor ^ is calculated with loop particles on mass shell: 

{k - q),^{k - q)u 



(4.2) 
(4.3) 

(4.4) 



J d(f)N^{kotr = qo, kotr - qo = 0) 

(4.5) 
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y^"" %-qy±2z{k^-qo)e 
-2^ + 4m - (1 - a){-fy + m 



{^-m-k')q' 



{k-qY±2i{ko-qo)e, 
In calculating the term proportional to (1 — a), we have to use the trick 

d 



(4.6) 



((A; - qY±it)-^ = lim 



{{k~qY±it)-m') 



2\-l 
7 



(4.7) 



For 7^ 0, we can decompose the vector k as 



2 -q + -^q + kT = {q- j^q)- 

r r \q\ 



2g2 



+ j^q + kr, 



(4.8) 



where, in the heat-bath frame with the z axis oriented along the vector q, we have 

g=(go,0,0,|g1), g = (|g1, 0, 0, g,), gg = 0, g' = -g^ (4.9) 
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The transverse component of k, kx vanishes after integration over 0. 
Finally, we obtain 



ijl + m)^(Jj + m) = 2m{q + m + 2m^) 



(4.10) 



2g2 > » ■ 

Now we can study the limit 

If{qo) = lim /^^(g^m^go) 

= (;^ + m) ° ^ N^{ko tr,ko tr - qo)- (4.11) 

It is easy to find that If{qo) is finite provided that 7^ and Ai'f^(q'o,0) < 00. The last 
condition is easy to investigate using the limiting procedure: 

Nniqo, 0) = lim iVn(^o, - qo) = lim 2n^{ko - qo)Mqo) - n^iko)) 

+ lim ((ne(go) - ne{ko) - {e{qo)e{ko - qo){ne{qo) + ^e(fco) - 2ne{qo)ne{ko))) ■ (4.12) 

One should observe here that the integration limits imply that the limit ko — > qo is taken 
from below for > m^, and from above for < m^. The two limits lead to different 



values of NQ{qo,0). Only the first term in (^.12|) can give rise to problems. We rewrite it as 



lim^^^o {^kon^{ko) ^ As relation ( [4.11|) should be vahd at any go, we can integrate 



over qo to find that the photon distribution should not grow more rapidly than \ko\ ^ as /Cq 
approaches zero, while the derivative of the electron distribution ne{qo) should be finite at any 
qo- 

lim koU^iko) < 00, (4.13) 

feo^O 

I^^Koo. (4.14) 
oqo 

Under the very reasonable conditions ( [4.13|) and ( [4.14| ) the electron propagator is free from 
pinches. 

It is interesting to observe the discontinuity of ijf (g^, m^, go) at the point g^ = m? . This 
feature will be repeated in massless QCD. 

It is worth observing that If{qo) is gauge independent, at least within the class of covariant 
gauges. 

4.2 Pinch Singularities in the Photon Propagator 

To consider the pinching singularity appearing in a single self-energy insertion approximation 
to the photon propagator, we have to make the substitutions mo = m = nis, smh."^ 9 ^{ko) — *■ 
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—ne{ko), smli^ 6s{ko — qo) — >■ —rieiko — qo), and h{ko) = —e{ko){l + 2n^{ko)). There are two 
thresholds, but only g^^^ = and the domain where < are relevant to a massless photon. 
The integration limits are given by the same expression (|4.2| ), but now we have to integrate 
over the domain (—00, ko 1) U(^o 2, +00). As q^ —0, we find {ko 1 —00) and {ko 2 — * +00). 
The integration domain is still infinite but is shifted toward ±00 where one expects that the 
particle distribution vanishes: 

K^u{q , g ) = U^p - (1 - a) — ] — \g^^-{l-a) ' 



g2 - 2iqoe ) q^ \ q^ + 2iqoe J 

+ f d^l^^niko, fco - qo) Lp -{I -a)- ^^^^ 



167r^\q\q^ \J-oo / \k\ -J ' \ q^ - 2iqoe J 

F'Hqo, fco, \klqk, .-) - (1 - «)^^^|— 1 ■ (4-15) 

\ q^ + 2iqoe J 

To calculate F^^" for the e — e loop, we parametrize the loop momentum k by introducing an 
intermediary variable / perpendicular to q. m is the mass of loop particles: 

k = aq + l, q.l = o, ^ = {k - qf = m\ = - a = + ' ' l)\ (4_ig) 

Iq 

At the end of the calculation we eliminate / in favor of k. After all possible singular denominators 
are canceled, one can set a = 1/2. 

F^^ = -Tr{jt + m)Y{j^ - k + 



+ C((4fco(A:o - qo) - ^m' - q')A^%q) + {-8{ko - |)' + 2g-«)i?^'^(g))). (4.17) 



Using relation (|7.21|) we obtain 



2 If f'^oi r°°\ kodko 



, <lo) = + J J I d^Nnik,, k, - g.) 



2^2 



4m q^ 



q^ 



2 

+ ^-({Ak,{k, - g„) - 4m' - g^) ^(g) + (-8(A;„ - |)2 + 2g-^)i?,,(g) ) ). (4.18) 



In the integration over kg the terms proportional to {k^q'^)" dominate and lim^2^Q \K^^{q'', qo) \ < 
00 if 
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Here N^{ko,ko- qo) is given by 

NQ,{ko, ko - Qo) = -2ne{ko - qo){-n^{qo) - n^iK)) 

-n^{qo) - rieiko) - e{qo)e{ko - qo){-n^{qo) + n^{ko) + 2n^{qo)ne{ko)) . (4.20) 

Assuming that the distributions obey the inverse-power law at large energies rieiko) oc |/co|~'^'' 
and ng{ko) oc Ifcol"*^*^, we find that the terms linear in densities dominate. Thus, for n = 0, 1, 
one finds 

q^ \J-oo Jko 2 / \k\ 



It follows that (|4.19|) is finite (in fact, it vanishes) if (5e, 5e > 3. Similar analysis for electron 



propagator at < (thus outside of our analysis of pinch singularities) leads to 5^ > 3. This 
is exactly the condition 

) < oo. (4.22) 

Thus the pinching singularity is canceled in the photon propagator under the condition that the 
electron and positron distributions should be such that the total number of particles is finite. 

Also, in the photon propagator, the quantity \img2 ^QK^^{q'^,qo) does not depend on the 
gauge parameter. 

Expression ([4.21|) is not valid for m = 0. 



5 Pinch Singularities in Massless QCD 

In this section we consider the case of massless QCD. Pinching singularities, related to massive 
quarks, are eliminated by the methods used in the preceding section. 

In self-energy insertions related to gluon, quark, and ghost propagators, the masses in the 
loop as well as the masses of the propagated particles are zero. Thus, the methods of the 
preceding section do not produce the expected result. Attention is turned to the spin degrees 
of freedom, i.e., to the function F of the integrand in ( |3.4| ) to ( p.9| ). In the calculation of F it 
has been anticipated that the loop particles have to be on mass shell. In this case, F provides 
an extra q^ factor in all the cases considered, in which not all particles are scalars. This 
factor suffices for the elimination of pinching singularities. 

The integration limits are now 

ko 1,2 = ^ (qo T W ■ (5.1) 

The difference \k\2 — \k\i is finite and there is no threshold effect. 

It is worth observing that for > 0, we have to integrate between koi and ko2, whereas 
for < 0, the integration domain is {— oo, ko i)[j{ko 2, +oo). This leads to two limits, 
lim^2^_i_o A'(g^, go) = K±{qo), in all cases of massless QCD. 

By inspection of the final results ( p.4[) , (|5.5|) , and (|5.6|) , we find that the case g^ < requires 
integrability of the function koNQ{ko,ko — qo) leading to the condition (|4.22| ) on the quark, 
gluon, and ghost distribution functions. 



By using ( [4.16[ ), we again introduce the intermediary variable I perpendicular to g; now we 



have to set m = 0. 
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5.1 Self-Energy Insertions Contributing to the Gluon Propagator 

The function K^j_y{q^ , qo) related to the gluon propagator is the sum 

'h^iKq^g^rnassless,fj.u{Q i 1o) ~^ ghgh,^v{,Q. ) Q.o) ~\~ K gg^^yi^q , qo) , (5-2) 

where the terms in the sum are defined as 

K^uiq'^, Qo) = {g^,p - (1 - a)DR^p{q)) — {g^^ - (1 - a)DAauiq))- (5.3) 

Pinching singularities, related to massive quarks, are eliminated by the methods used in the 
preceding section. The tensor F related to the massless quark- ant iquark contribution to the 
gluon self-energy is 



^ , - Qo) - q^)A^^{q) + (-8(A;, - ^-^)' + 2<f)B>^^{q) ] + O^^ {kT)] .{^A) 



\q^ \ 2 



As F^jy contains only A and B projectors, relation (|7.21|) guarantees that the result does not 
depend on the gauge parameter. 

Relation (|5^) contains only terms proportional to g^, and limg2^Q Kfj,iy{q'^ , qo) is finite. 

For the ghost-ghost contribution to the gluon self-energy, the tensor F is given by 



-<5a.iv4( ^^°^^" g^°^^^' A"'^(g) - (ko - ^)'B^'^iq) - ^D^^{q) + 0'^^(m)(5.5) 



(ff^ko{ko-qo) + q^^^y,^, qo^2T^^u^^^ f 

r 

The antisymmetric part vanishes after integration, so we have left it out from the final result 



in (U). 



The tensor F for the gluon-gluon contribution to the gluon self-energy is 

K = ^ (a'^-iq + ky - g''^{2k - q)^ + g^>^{k - 2q)^) 

( /I ^ {k - q)a{k - q)a' \ 

gaa' - (1 - a' 



{k-q)^±2t{ko-qo)eJ 
(^-'(g + k)-' - g-'^'{2k - g)^ + /^(fc - 2g)'^') (^g,,, - (1 - < 



A;2 ± 2ikoe 



2 V 2 

-1,2 4 

1 — a f q 27 7 2 /7 1 \ q 

-rqfjiqv + 5g Ipl^ + q [l^qy + l^q^) + —gf,^ 
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— CL / 5' r 27 7 2 /7 7 N y 



+ (-10(fc, - |)2 + 4g-^)i?'^^(g) - ^^^'^(g) 
-(1 - a) ( -A^^ - B^"" - ^C"" 

2 ^2 



+ (1 - a)2 + 2^5'^'^ - 2^C"^^ - 2D'^^) + O^'^(A^r)) . (5.6) 

Expressions (|5l^ ), ( ^.5|) , and (|5.6|) have been obtained by substitution of (|7.17|) , (|7.19|) , and 
( [r.2(]| ) and, finally, by eliminating the variable I in favor of k. The tensor 0^'^{kj') is linear in 
k, thus it vanishes after integration over 0. 

We note here that, in the Feynman gauge (a = 1), the operator C is absent from the gluon 
self-energy! Consequently, the relation originating from Slavnov- Taylor identities (proved in 
2^ , |25| , |26[] for equilibrium densities), 11^ = (g^ — IIl)!!/), is fulfilled at a = only if = 0. 
Thus the contributions to ttd from the ghost-ghost and gluon-gluon self-energies mutually 
cancel, imposing restrictions on the densities related to unphysical degrees of freedom. As it 
does not interfere with the cancellation of pinches, the problem of unphysical degrees of freedom 
will be discussed elsewhere. 

Finally, we need ( [7.21|) in all three cases. 

Expressions ( ^^) , (|5.5| ), and ( ^.6|) for the ghost-ghost, quark- antiquark, and gluon-gluon 
contributions to the gluon self-energy contain only terms proportional to g^. The function 
K^y{q^,qo) approaches the finite value K^y{±,qo). 

Thus we have shown that the single self-energy contribution to the gluon propagator is free 
from pinching under the condition (|4.22|) . 

5.2 Quark-Gluon Self-Energy Contribution to the Quark Propaga- 
tor 

The K spinor for the quark-gluon contribution to the massless quark propagator is defined as 

H{q\qo)=li%ll- (5.7) 
In the self-energy of a massless quark coupled to a gluon the spin factor F is given by 

w A -1 ( , (fc-g)^(fc-g). ^1 , 

= ^"^^iVT -^^""hk- qy ± 2.(fc. - q^)e) ^ 
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'^afc -2^ - (1 - a)(-^ 



For our further discussion, we need the product 

= Sa,^^^' (-^ - 1^ + + - ^(-^ + 1^)) , (5.9) 

which contains the damping factor q^. The term ^j, will be integrated out. 
By inserting ( p.9|) into (|5.7| ), we obtain ( 2.37|) free from pinches. 
To calculate If{qo), we need the limit 

hm — f- = ^afe^TTT (5.10) 



From ( |5.1CI| ) we conclude that If{qo) does not depend on the gauge parameter. 



Omitting details, we observe that pinching is absent from the quark propagator, also in the 
Coulomb gauge. 

5.3 Ghost-Gluon Self-Energy Contribution to the Ghost Propaga- 
tor 



The K factor is defined as 



K{q\q,) = %. (5.11) 



The F factor for the ghost-gluon contribution to the ghost self-energy is 

(^M - - 



5abN,k^^q-\g^,-{l~a) 



{k-qy±2i{ko-qo)e, 



^ 5abN,'^. (5.12) 

The factor q^ ensures the absence of pinch singularity and a well-defined perturbative result. 

5.4 Scalar-Photon Self-Energy Contribution to the Scalar Propa- 
gator 

Although the massless scalar boson interacting with a photon is not part of massless QCD , it 
is treated using the same methods. 
The K factor is defined as 



K{q\q,) = -^. (5.13) 



q^ 
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The F factor for the massless scalar-photon contribution to the scalar self-energy, 
F,^ = {q + kriq + ky U^^-(l-a)- 



{k-q)^{k-q)^ 



{k - qf ± 2i{ko - qo)e, 



g 2 — (1 — a 



{k-qY±2i{ko-qo)ej 
2q\ (5.14) 



clearly exhibits the q^ damping factor! 



6 Conclusion 



Studying the out of equilibrium Schwinger-Dyson equation, we have found that ill-defined pinch- 
like expressions appear exclusively in the Keldysh component {Qk) of the resummed propagator 
( p.32| ), or in the single self-energy insertion approximation to it (|2.35|) . This component does 



not vanish only in the expressions with the Keldysh component ( |2.21| ) {Vt or Vt for the single 
self-energy approximation) of the self-energy matrix. This then requires that loop particles be 
on mass shell. This is the crucial point to eliminate pinch singularities. 

We have identified two basic mechanisms for the elimination of pinching: the threshold and 
the spin effects. 

For a massive electron and a massless photon (or quark and gluon) it is the threshold effect 
in the phase space integration that produces, respectively, the critical q^ — im? or q^ damping 
factors. 

In the case of a massless quark, ghost, and gluon, this mechanism fails, but the spinor/tensor 
structure of the self-energy provides an extra q^ damping factor. 

We have found that, in QED, the pinching singularities appearing in the single self-energy 
insertion approximation to the electron and the photon propagators are absent under very 
reasonable conditions: the distribution function should be finite, exceptionally the photon 
distribution is allowed to diverge as k~^ as ko 0; the derivative of the electron distribution 
should be finite; the total density of electrons should be finite. 

For QCD, identical conditions are imposed on the distribution of massive quarks and the 
distribution of gluons; the distributions of massless quarks and ghosts (observe here that in 
the covariant gauge, the ghost distribution is not required to be identically zero) should be 
integrable functions; they are limited by the finiteness of the total density. 

In the preceding sections we have shown that all pinch-like expressions appearing in QED 
and QCD (with massless and massive quarks!) at the single self-energy insertion level do 
transform into well-defined expressions. Many other theories behave in such a way. However, 
there are important exceptions: all theories in which lowest-order processes are kinematically 
allowed do not acquire well-defined expressions at this level. These are electroweak interactions, 
processes involving Higgs and two light particles, a p meson and two vr mesons, Z, W, and other 
heavy particles decaying into a pair of light particles, etc. The second important exception is 
massless g'^(f)^ theory. This theory, in contrast to massless QCD, contains no spin factors to 
provide (|1.6|). In these cases, one has to resort to the resummed Schwinger-Dyson series. 
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The main result of the present paper is the cancellation of pinching singularities at the 
single self-energy insertion level in QED- and QCD-like theories. This, together with the 
reporte cancellation of coUinear singularities, allows the extraction of useful physical 

information contained in the imaginary parts of the two-loop diagrams. This is not the case 
with three-loop diagrams, because some of them contain double self-energy insertions. In this 
case, one again has to resort to the sophistication of resummed propagators. 
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7 Appendix 

We start [0 by defining a heat-batli four- velocity U^, normalized to unity, and define the or- 
thogonal projector 

\u = g^.u - U^U,. (7.1) 
We further define spacelike vectors in the heat-bath frame: 

= A^,g", /t^/t'^ = K^ = -f. (7.2) 

There are four independent symmetric tensors (we distinguish retarded from advanced ten- 
sors by the usual modification of the ie prescription) A, B, and D (which are mutually orthog- 
onal projectors), and C: 

A,M) = \. - (7.3) 

r (n)-t^f^(^^riri "^m^'- , ^o + ^ \ ,^ 



+ 2iqoe 



+ 2 

In addition to the known multiplication!!' properties 

A{q)A{q) = A{q), BR^A{q)BRAq) = BrAq), (7-7) 

CR,AiQ)CR,AiQ) = -(BrAq) + DrAq))^ DRAq)DRAq) = DrAq), (7.8) 
A{q)B{q) = B{q)A{q) = 0, A{q)C{q) = C{q)A{q) = 0, 

Aiq)Diq) = D{q)A{q) = 0, i?(g)D(g) = Diq)Biq) = 0, (7.9) 

{BRAq)CRAq)),,^ = {CRAq)DRAq)),u = , "^f:. , (7.10) 



{CRAq)BRA<i)),u = {DRAq)CRAq)),u = , "tt , (7.ii) 

we need mixed products 



BRAq)BAAq) = liBniq) + BA{q)), (7.12) 
CRAq)CAAq) = -^iBRiq) + BAiq) + DR{q) + DA{q)\ (7.13) 
DRAq)DAAq) = liDRiq) + D^(g)), (7.14) 



[BRAq)CAAq)),^ = {CRAq)DAAq)),^ = ^(^4%^ + q'-l^qj ' ^^'^^^ 
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{CnA<l)BAA<l)),^ = iDnAl)CAAl)),u = T^i '',. + , '"Z. J - (7-16) 



By calculating the traces of the tensors IH'^ , q^^l'^ + Pg*^, and q^q'^ with projectors, we find 



qn" + Pg^ = -^C^''(g) + O^^'ikr), (7.18) 



= m^^A^-(g) + 1_ ( ^A^'iq) - liB^^iq) I + O^^ikr), (7.17) 



k1 



qi'q'^ = q^D'^'iq), (7.19) 

= {A^^'iq) + B^'^iq) + D^'^^q)). (7.20) 

The tensor O^^ijir) is linear in fc, and vanishes after integration over (j). One should observe 
that ( |7.17| ) to ( [7. 201 ) are valid for an arbitrary (but the same for B and D) R/A prescription, 
so we do not indicate it. 

Using the multiplication rules one obtains 

{g^p - (1 - a)Dn,p{q)){UAP^ + hB^^ + f.C^'^ + f^D^'^ig,, - (1 - a)DA.u{q)) 
= I (/a< + hB^R - (1 - «)/cCr + (1 - aff^D^n) 

+l(faA7 + hB^/ - (1 - a)/,Cr + (1 - aff.D^/). (7.21) 



It is important to observe that, owing to the properties of the mixed products ( [7.12| ) to (|7.16| ), 
the R/A assignment of F does not infiuence the final result! 
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